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A NOTE ON GENERALIZED MALLIAVIN CALCULUS 

S. V. LOTOTSKY, B. L. ROZOVSKII, AND D. SELESI 



Abstract. The Malliavin derivative, divergence operator, and the Ornstein- 
Uhlenbeck operator are extended from the traditional Gaussian setting to gener- 
ahzed processes from the higher-order chaos spaces. 



1. Introduction 



Stochastic integration started with the construction of integrals with respect to the 
Wiener process [^ and then extended to a much larger class of processes [17] . The 
Wiener process was also in the core of the early development of Malliavin calculus 
[To], but generalizations so far have not been nearly as sweeping as in the theory of 
stochastic integration. Currently, the driving random source in Malliavin calculus is 
an isonormal Gaussian process on a separable Hilbert space [HI [13]. This process, 
similar to the standard Wiener process, is in effect a linear combination of a countable 
Q\ ■ collection ^ := {^j}j>i of independent standard Gaussian random variables. 

^ : 

lO ■ A natural question to ask is whether one can extend Malliavin calculus to nonlinear 

^ ! functionals of isonormal Gaussian process as the driving random source, and still enjoy 

r^ I all the benefits of the Gaussian setting. Natural candidates for this role are elements 

Q ' of the Hilbert space of square integrable functionals of the isonormal Gaussian process. 

This space is often referred to as Wiener Chaos space. 

In this paper we extend the main operators of Malliavin calculus to the space of 
/\ I generalized random elements X]|a|<oo /a^"' where {^q, |q:| < oo} is the Cameron- 

c^ ■ Martin basis in the Wiener Chaos space, a is a multiindex and fa belong to a certain 

Hilbert space X (for detail see Section 2). To cover some emerging applications, we 
allow formal linear combinations with infinite variance, that is X]|a|<oo ll/allx ~ ^^• 

Looking for solutions that are generalized random elements is quite reasonable: after 
all, the Gaussian white noise that often drives the equation of interest is itself a 
generalized random element. Our interest in this subject was prompted by some 
recent and not so recent developments in the stochastic partial differential equations 
(SPDEs). These developments indicate that large classes of solutions of linear and 
nonlinear SPDEs driven by Gaussian sources are generalized random elements. 
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One example is the heat equation driven by multiphcative space-time white noise 
W{t,x) with dimension of x two or higher [1] 



2 



ut = Au + uW, u{0, x) = e-l^l^ (1.1 

Examples in one space dimension also exist: a stochastic parabolic equation violating 
the parabolicity condition [7]: 

du = Uxxdt + aUxdw(t), cr^ > 2, u{0,x) = e~^ , (1.2) 

or a stochastic parabolic equation of full second order [8]: 

2 

du = Uxxdt + Uxxdw{t), u{0,x) = e~^ . (1.3) 

In all three examples, E||M(t, ■)|||^ = J2\cx\<oo \Wa(t,-)\\'x = ^^ ^^ all t > and all 
typical function spaces X, such as Sobolev spaces. 

As a different example, consider equation 

- (a (x) Uxix))^ = fix), X e (0, 1) , M (0) = M (1) = 0, (1.4) 

with a{x) = d{x) + e(x), where d{x) is non-random and e{x) = J2k>i'^k{x)^k is a 
Gaussian noise term; J2k>i ^^Px ^fc(^) < ^^- Recently, this equation was investigated 
in the context of uncertainty quantification for mathematical and computational mod- 
els [19]. As problem (II. 4p is ill posed, one could modify it as follows: 

- (a (x) Vx{x))^ + {6,^x) {vx (x)))^ = /(x), , . 

xG(0,l),t;(0) = t;(l) = 0, ^ ' 

where Se(^x) stands for Malliavin divergence operator (Skorokhod integral) with respect 
to Gaussian noise e(x). In contrast to (II. 4p . equation (11.51) is well posed and uniquely 
solvable, and similar to equations (ll.ll) - (ll.3p . E ||f H^^ = oo for all traditional spaces X 
of functions on (0, 1). Technically, the above modification of problem (11.41) amounts 
to replacement of products of random elements by stochastic convolutions, such as 
Wick products[3llll[20]. In the literature on quantum physics, procedures of this type 
are often called stochastic quantization [2], [18]. Equations subjected to the stochastic 
quantization procedure are usually referred to as quantized. 

We remark that the replacement of equation (ll.4p by equation (ll.Sp also mimics the 
idea of Ito [3] of replacing the singular equation 

x{t) = a (x (t)) + cr (x (t)) w{t) 

by the well posed stochastic differential equation 

X (t) = xo + / a{x{s))ds+ / a {x (s)) dw{s). 
Jo Jo 

Because there is no natural filtration associated with elliptic equation (11. Sp . the Ito 
integral has to be replaced by the Skorokhod integral. 

Equation (II. 5p differs from (II. 4p quite drastically. While both equations are stochastic 
perturbations to the solution of the deterministic equation 

- (a (x) Vx{x))^ = fix), X G (0, 1) , t; (0) = ^ (1) , = (1.6) 
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only the solution to the quantized equation fll.Sp is an unbiased perturbation of the 
solution of equation fll.6p in that Ef (x) = v (x) is a solution of equation (ll.6p : the 
solution of equation (11 ■4p . even if existed, would not enjoy this property. 

Two other examples of stochastic quantization are currently under investigation: ran- 
domly forced Burgers equation [6] and Navier-Stokes equation [12]. Let us consider 
Burgers equation 

Ut = u^x + uu^ + e"^ ^, t > 0, X eM, (1.7) 

with a deterministic initial condition, where ^ is a standard Gaussian random variable. 
The stochastic quantization of this equation is 

vt = Vxx + Sv (vx) + e""^'^, (1.8) 

where 6^ (vx) is the Malliavin divergence operator (Skorokhod integral) of Vx with 
respect to the solution v of (II. 8p (and v is not a Gaussian process). It is shown 
[6] that V {t,x) := Kv {t,x), with a suitable interpretation of the expectation, is the 
solution of the deterministic Burgers 

Vt (t, x) = v^^ (t, x) +v (t, x) v^ (t, x) . 

Thus, as in the linear example (II. 4p . the quantized version of stochastic Burgers 
equation (II. 8p is an unbiased perturbation of the deterministic Burgers equation with 
the same initial condition. The standard stochastic Burgers equation does not have 
this convenient property. Similar effect holds for quantized Navier-Stokes equation 

m. 



We emphasize that, in all examples we have discussed, the variance of a generalized 
random element u is infinite and is given by the diverging sum Y2a II'^qIIx ~ ^^• 
However, the rate of divergence differs substantially from case to case. To study 
this rate of divergence, we introduce a rescaling, or weighting, operator 91 defined by 
9^a = f'ct^ai where weights r^ are positive numbers selected in such a way that the 
weighted sum ^^ r^ ||tia||j(^ becomes finite. Of course, the choice of r^ is not unique 
and depends on the specifics of the problem, for example on the type of the stochastic 
PDE in question. A special case of this rescaling procedure originates in quantum 
physics and is related to second quantization 



Quantum physics has brought about a number of important precursors to Malliavin 
calculus. For example, creation and annihilation operators correspond to Malliavin 
divergence and derivative operators, respectively, with respect to a single Gaussian 
random variable. The original definition of Wick product [20] is not related to the 
Malliavin divergence operator or Skorokhod integral but remarkably these notions 
coincide in some situations. In fact, standard Wick product could be interpreted as 
Skorokhod integral with respect to square integrable processes generated by Gauss- 
ian white noise, while the classic Malliavin divergence operator integrates only with 
respect to isonormal Gaussian process. In Section [31 we demonstrate that Malliavin 
divergence operator can be extended to the setting where both the integrand and the 
integrator are generalized random elements in a Hilbert space, although we did not 
try to extend Wick product in a similar way. 
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In this paper, we restrict ourselves to the basic study of the three main operators 
in the Malhavin calculus: the derivative operator D, the divergence operator S and 
the Ornstein-Uhlenbeck operator C = S o D. We present constructions of D„(w), 
^m(/); Ei'iid Cu{v) when u,v, f are Hilbert space-valued generalized random elements. 
Section [2] reviews the main constructions of the Malliavin calculus in the form suitable 
for generalizations. Section |3] presents the definitions of the Malliavin derivative, 
Skorokhod integral, and Ornstein-Uhlenbeck operator in the most general setting of 
weighted chaos spaces. Section H] presents a more detailed analysis of the operators 
on some special classes of spaces. 

To illustrate some of the main results, let us consider the one-dimensional setting. 
Let ^ be a standard normal random variable and define 



where H„ is nth Hermite polynomial. If / is a square-integrable functional of ^, then 

n>0 

The space of square-integrable functionals of ^ can thus be identified with £2'- 

{/n, n>0:^/2<oo}. 

n>0 

We define a generalized random functional / of ,^ as a collection of numbers {/„, n> 
0} without any restrictions on /„ and a formal representation 

fc>0 

Let M, /, V be generalized functionals of ^ and let p, g, r be positive real numbers such 

that 

1 1 _ 1 

p q r ' 

We show in the paper that if 

n>0 n>0 



then Y)u{v) is a generalized functional of ^ such that 



and 



Similarly, if 



fc=o ^ • • / 



n>l ^ 



n>0 n>0 
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then du{f) is a generalized functional of ^ such that 6u{f) = uo f, where operator o 
stands for the Wick product, 






and 

J2r-{SM)l <oc. 

n>l 

Finally, if p, q, r are positive real numbers such that 

r pj \ pj 
(for example, p = 1, q = '2, r = 1/2) and 

n>0 n>0 

then Cu{v) is a generalized functional of ^ such that 

n oo 
k=0 m=0 

and 

n>0 

2. Review of the traditional Malliavin Calculus 

The starting point in the development of Malliavin calculus is the isonormal Gaussian 
process (also known as Gaussian white noise) W: a Gaussian system {W{u), u G lA} 
indexed by a separable Hilbert space lA and such that ¥W{u) = 0, K(\V{u) W{v)) = 
(m, v)u- The objective of this section is to outhne a different but equivalent construc- 
tion. 

Let F := (fi, J-", P) be a probability space, where J-" is the cr-algebra generated by 
a collection of independent standard Gaussian random variables {^i}j>i. Given a 
real separable Hilbert space X, we denote by L2{F]X) the Hilbert space of square- 
integrable J-"- measurable X- valued random elements /. When X = M, we often write 
L2(F) instead of L2(F;M). Finally, we fix a real separable Hilbert space U with an 
orthonormal basis il = {u^, k > 1}. 

Definition 2.1. A Gaussian white noise W onU is a formal series 

W = J2^kUk. (2.1) 

fe>i 
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Given an isonormal Gaussian process W and an orthonormal basis it in U, represen- 
tation (12. ip follows with ^k = W{uk). Conversely, (12 .ip defines an isonormal Gaussian 
process on U by 

fc>i 
To proceed, we need to review several definitions related to multi-indices. Let JT" 
be the collection of mult i- indices ex = (ai, 02, . . .) such that a^ G {0, 1, 2, . . .} and 
J2k>i '^k < 00. For OL, (3 E J , we define 

q; + /3 = (ai + /3i,a2 + /32,---)5 l«l = X^"fc, Q:! = ]^afc!. 

fc>i fc>i 

By definition, a > if |q;| > and /3 < ck if Pk < otk for all /c > 1. If /3 < ck, then 

ex - 13 = {ai - I3i,a2 - P2, ■ ■ ■)■ 
Similar to the convention for the usual binomial coefficients, 

( a' 

'ol\ 7 ^77^, if /3 < a. 



(;)^ («-/3)i.r 



0, otherwise. 

We use the following notation for the special mult i- indices: 

(1) (0) is the multi-index with all zero entries: (0)^. = for all /c; 

(2) £{i) is the multi-index of length 1 and with the single non-zero entry at position 
i: i.e. £{i)k = 1 ii k = i and s{i)k = if /c 7^ i. We also use convention 
s(0) = (0). 

Given a sequence of positive numbers q = (gi, g2, • • •) and a real number £, we define 
the sequence q^"", ex ^ J ,hy 

k 

In particular, 

(2N)^° = JJ(2A;)^°^ 
fc>i 
Next, we recall the construction of an orthonormal basis in L2(F;X). Define the 
collection of random variables 

H = {^c., exej] 
as follows: 



««=n 



where H„ is the Hermite polynomial of order n: 

Sometimes it is convenient to work with unnormalized basis elements Hq,, defined by 

H„ = y/^\U- (2.2) 
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Theorem 2.2 (Cameron-Martin [T]). The set S is an orthonormal basis in L2(F; X): 

II?; I|2 



if V E L2(F;X) and Va = "^{v ^a) , then v = J^aej^'^^^ ^''^^ ^^ll^llx = S-^ -^ ll^< "^ 



If the space U is finite-dimensional, then the multi-indices are restricted to the the 
set 

Jn = {oi e J : ak = Q, k > n}. 

The three main operators of the Malliavin calculus are 

(1) The (Malliavin) derivative D^^; 

(2) The divergence operator (5^, also known as the Skorokhod integral; 

(3) The Ornstein-Uhlenbeck operator C-^ = 6^rT)^r■ 

For reader's convenience, we summarize the main properties of D^ and ^^r', all the 
details are in [13], Chapter 1]. 



(1) D^ is a closed unbounded linear operator from L2(F; X) to L2{¥; X ®U)\ the 
domain of Y)y^ is denoted by ro^'^(F; X); 

(2) \i V = F[W{hi),...,W{hn)) for a polynomial F = F{xi,...,Xn) and 
hi, . . . ,hn E X, then 

fc=i ^'^ 

(3) (5yj/ is the adjoint of D^^ and is a closed unbounded linear operator from 

L2(F; X ® W) to L2(F; X) such that 

EMvK(/))=E(/,D^(y,))^ (2.4) 

for all ip G D^'2(F; M) and / G 3^'^{¥; X ® W). Equivalently, 

(^'><5h/(/))z,,(F;X) = (/'Dh/M),.,(F;x®w) (2-5) 

for all w G Di'2(F; X) and / G ©^'^(F; X ® W). 

We will need representations of the operators D^y, dy^, and £^y in the basis S. 

Theorem 2.3. (1) If v G i:2(F;X) and 

y^ |q:| IK'alli- < oo, (2.6) 

then Dyyiv) G L2(F; X ®U) and 

^wi'") = 5Z X] ^/a^ia-e{k) Vc ® Ufc. (2.7) 

aej fc>l 

a6j',fe>l 

anc? 

5^ |a|||/fc,,||^<oo, (2.8) 

aej,k>l 
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then Syy{f) e L2{¥;X) and 

^wU) = Y. ( Ev^^'^.^-W ) ^"- (2.9) 

aej \k>l J 

(3) Ifv eL2{¥;X) and 

J2H^\K\\x<oo, (2.10) 

then C^/{v) G L2{¥;X) and 

'CvyH = I]l"K'c.eo. (2.11) 

Proof. Linearity of the operators implies that, in each case, it is enough to find the 
image of ^a. 

(1) Using (12. 3 p and properties of the Hermite polynomials, for every ex. & J', 

Diy({a) = Y VokU-eik) Ufc. (2.12) 

fc>l 
Orthonormahty of {C,a, ct G J^} and {ufc. A; > 1} then implies 

JE||D^^(v)||^jgj^ = ^ A/afc /3n IE(^a-£(fe)^/3-6{n)) (Vo,, V,3)x (Ufc, U„)w 
ct,l3,k,n 

Ell l|2 Y^ I I II ||2 

ci,k a 

(2) By (1212|) and (EH), for every ^„ G H, /i G X, and Ufc G il. 



5^(/i(g)Ufc^a) = h^Jak + l^Q+Efc. (2.13) 

(3) By (ITT^ and (ITT^ . for every ^«, 

-CvKiec) = <5^(D^^(e«)) = |a|^„. (2.14) 

D 

Remark 2.4. Here is an important technical difference between the derivative and 
the divergence operators: 

• For the operator D^, 

{T>^{v))^ = ^ V«fc + 1 Vc+e{k) ® Uk] (2.15) 

fc>i 

in general, the sum on the right-hand side contains infinitely many terms and 
will diverge without additional conditions on v, such as (12. 6p . 

• For the operator dyy, 

fc>i 
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the sum on the right-hand side always contains finitely many terms, because 
only finitely many of ak are not equal to zero. Thus, for fixed a, (<5^(/)) is 
defined without any additional conditions on /. 



3. Generalizations to weighted chaos spaces 



Recall that W, as defined by (12. ip . is not a W- valued random element, but a generalized 
random element on U: 

W(/i) = ^(/i,Ufc)w6, (3.1) 

k>l 

where the series on the right-hand side converges with probability one for every h ElA. 
The objective of this section is to find similar interpretations of the series in (12. 7p . 
(EHD, and (EUD if the corresponding conditions ([2SD, (E3D, fl^^ fail. Along the 
way, it also becomes natural to allow other generalized random elements to replace 
W. 

We start with the construction of weighted chaos spaces. Let 7^ be a bounded linear 
operator on L-^iW) defined by TZ^a = fcxia for every a G JT", where the weights 
{tq,, a G J} are positive numbers. 

Given a Hilbert space X, we extend TZ to an operator on L2(F; X) by defining TZf as 
the unique element of Iv2(F; X) so that, for all g G Iv2(F; X), 

E(7^/,^)x = 5^r«E((/,^)xeo)• 
Denote by TZL2{¥\ X) the closure of L2(F; X) with respect to the norm 

ll/ll7^L2(F;X) •= II'^/IIl2(F;X)- 

In what follows, we will identify the operator TZ with the corresponding collection 

(ra, OL E J). Note that if m G 7?.iL2(F; X) and v G 7?.2L2(F; X), then both u and v 
belong to 7?.L2(F;X), where r^ = min(ri_Q,r2,a)- As usual, the argument X will be 
omitted if X = M. 

Important particular cases of 7^L2(F;X) are 

(1) The sequence spaces L2.q(F;X), corresponding to the weights 



where q = {qu, A; > 1} is a sequence of positive numbers; see O [3 [H]. Given 
a real number p, one can also consider the spaces 

L|^(F;X) = L2,q.(F;X), (3.2) 

where q^ = {g^, /c > 1} . In particular, L\^ = L2,q;L^q = L2,i/q- Under the 
additional assumption gfc > 1 we have, similar to the usual Sobolev spaces, 

LU¥;X)cLU¥;X),p>r. 
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(2) The Kondratiev spaces {S)p/{X), corresponding to the weights 

r« = (a!)^/2^2N)^", p G [-1, 1], i G M, (3.3) 

see [1]. 

There is a natural duahty between L2,q(X) and L^q(X): 

{u,v)^=^{uc„v„)x] (3.4) 

there is a natural duality between {S)p/{X) and (iS)_p^_^(X): 

(m, t;)p,£ = ^ (mq, t^«)x- (3.5) 

Both (13.41) and (13.51) extend the notion of K{u,v)x to generalized X-valued random 
elements. 

Taking projective and injective limits of weighted spaces leads to constructions similar 
to the Schwartz spaces iS(]R'^) and iS'(]R'^). Of special interest are 

(1) The power sequence spaces 

L+^(F;X) = f|L|^(F;X), L^^^{¥;X) = [j Ll^{¥;X), (3.6) 



where q = {qk, k > 1} is a sequence with g^ > 1 (see (13.21) ). 
(2) The spaces S''{X) and 5_p(X), < p < 1 of Kondratiev test functions 
and distributions: 

S^{X) = []{S)pAX), 5_p(X) = \J{S).pAX). (3.7) 



In this regard we mention that, in the traditional white noise setting, X = M*^, 
p = corresponds to the Hida spaces, and the term Kondratiev spaces is 
usually reserved for iS"^(M'^) and 5_i(M'^). 

If the space U is finite-dimensional, then the sequence q can be taken finite, with as 
many elements as the dimension of U. In this case, certain Kondratiev spaces are 
bigger than any sequence space. 

Proposition 3.1. IfU is finite- dimensional, then 

L2,q(X) C (5)_p,_,(X) (3.8) 

for every p > 0, i > p and every q. 

Proof. Let n be the dimension of U and r = min{gi, . . . , g„}. Define r = {r, . . . , r}. 
Then L2,q(X) C L2,t(^)- On the other hand, for all a e J, {2N)^cx\ > \cx\\, 

(r2|°l(a!)''(2N)2'')-i < (r2|"l(|a|!)0"' < C(r) 
and therefore L2,r{X) C (5)_p,_£(X). D 
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Analysis of the proof shows that, in general, an inclusion of the type f l3.8p is possible 
if and only if there is a uniform in a bound of the type (q^"(|Q:|!)'')^^ < C(2N)p"; 
the constants C and p can depend on the sequence q. If the space U is infinite- 
dimensional, then such a bound may exist for certain sequences q (such as q = N), 
and may fail to exist for other sequences (such as q = exp(N). Thus, both L2,q(X) 
and {S)-p^e{X) can be of interest in the study of stochastic differential equations. 

Definition 3.2. A generalized X-valued random element is an element of the 
set [JTZL2{¥] X) , with the union taken over all weight sequences TZ. 

To complete the discussion of weighted spaces, we need the following results about 
multi-indexed series. 

Proposition 3.3. Let r = {r^, k > 1} be a sequence of positive numbers. 
(V //'Efc>i^fe < oo, then 

aej ' \k>l / 

(2) IfYlk>i '"*: ^ '^^ '^'^^ ^k < ^ for all k, then, for every oc ^ J , 

/3ej \ /^ / \fc>i k J 

where \ — x is the sequence {1 — r^, A;>1}. /n particular, 

^(2N)^^"<oo (3.11) 

for all I > 1; cf jS Proposition 2.3.3]. 

(3) For every a E J', 

J2('^y = (^ + '^"^ (3.12) 

where 1 + t is the sequence {1 + r^, k > 1}. 
Proof. Note that 

\k>l / k>l n>l k>l fc>l 7i>l 

By assumption, limfc_^oo ^fc = 0, and therefore 

fc>i 

unless only finitely many of Uk are not equal to zero. Then both (13. 9p and (13.101) with 
CK = (0) follow. For general a, (I3.10p follows from 

Efn + k\ ^ 1 I I 1 

fc>0 ^ -^ ^ ' 
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which, in turn, follows by differentiating n times the equality J2k^'' ~ (1 — a^)~^- 
Recall that 

^rfc < cx), < rfc < 1 =^ 0<]^- < oo. 

k k ^'^ 

Equality fl3.12p follows from the usual binomial formula. D 

Corollary 3.4. (a) For every collection /q,, ol ^ J oj elements from X there exists 
a weight sequence r^, ex. ^ J such that 'Yla^^j ll/allx'^a < o*^- 

(^) If Qk > 1 (ind J2k>i ^/lk ^ °*^' ^^^''^ ^^^ space L^q(X) is nuclear, 
(c) The space S^{X) is nuclear for every p G [0, 1]. 

Proof, (a) In view of (13. lip , one can take, for example, 

(2N)-" 
'°"l+ll/.||x- 

(b) By (I3:T0|) . the embedding Lf+^(X) C L^^iX) is Hilbert-Schmidt for every peR. 

(c) Note that Efc>i(2^)"^ < oo. Therefore, by (l310D . the embedding (5)p,^+i(X) C 
{S)p,e{X) is Hilbert-Schmidt for every ieR. ' D 

To summarize, an element / of TZL2{¥;X) can be identified with a formal series 
Y^aej faia, whcrc /a e X and Y^aaj WfoWx^l < oo. Conversely, every formal series 
Tliaej fctia, /a G X, is a generalized X- valued random element. Using (12.21) . we get 
an alternative representation of generalized X-valued random elements: 

/=5^/«H«, (3.13) 

with f^ e X. 

The following definition extends the three operators of the Malliavin calculus to gen- 
eralized random elements. 

Definition 3.5. Let u = Ylaej'^a^a be a generalized U -valued random element, 
V = '^fj^^jVaC.a, 0- generalized X -valued random element, and f = '^^^j fa^a, o- 
generalized X ® U-valued random element. 

(1) The Malliavin derivative ofv with respect to u is the generalized X®U -valued 
random element 



provided the inner sum is well-defined. 

(2) The Skorokhod integral of f with respect to u is a generalized X -valued random 
element 

^«(/) = E (E J (^) (^'«-/3)"] U- (3.15) 
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(3) The Ornstein-Uhlenbeck operator with respect to u, when applied to v, is a 
generalized X -valued random element 

^^(^') = J2 [ J2 y [g) [ ^^1 ^'/3+7K.^"-/3)w j^o, (3.16) 

provided the inner sum is well-defined. 

The definitions imply that both D and 6 are bi-linear operators: 

^au+bv{u^) = aAu{w) + bAv{w), Au{av + bw) = aAu{v) + bAu{w), a, 6 G M, 

for all suitable u, v, w; A is either D or 6. The operation Cu{v) is linear in v for fixed 
u, but is not linear in u. The equality 

De,(ec.) = 

shows that, in general J^u{v) 7^ D^(u). The definition of the Skorokhod integral du{f) 
has a built-in non-symmetry between the integrator u and the integrand /: they have 
to belong to different spaces. This is necessary to keep the definition consistent with 
f l2.9p . Similar non-symmetry holds for the Ornstein-Uhlenbeck operator Cu{v). Still, 
we will see later that du{f) = ^f{u) if both / and u are real- valued. If D„(w) is defined, 
then Cu{v) = <5„(D„(t;)), but Cu{v) can exist even when Du{v) is not defined. 

Next, note that Su{f) is a well-defined generalized random element for all u and /, 
while definitions of D„(w) and Cu{f) require additional assumptions. Indeed, (^) = 
unless /3 < a, and therefore the inner sum on the right-hand side of (13.151) always 
contains finitely many non-zero terms. By the same reason, the inner sums on the 
right-hand sides of (13.141) and (I3.16P usually contain infinitely many non-zero terms 
and the convergence must be verified. This observation is an extension of Remark 
12.41 and we illustrate it on a concrete example. The example also shows that Cu{v) 
can be defined even when D„(w) is not. 

Example 3.6. Consider u = v = W. Then D„(f ) is not defined. Indeed, 

fufc, if a = e{k), k>l, 
Uoc = Vc. = \ . (3.17) 

I U, otherwise. 

Thus, {Du{v))^ = if |a| > 0, and 

fc>i 
which is not a convergent series. 

On the other hand, interpreting v as an M (g) W- valued generalized random element, 
we find 

,, , ,, fy2,c. = 2s(fc), k>l, 
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or, keeping in mind that \/2^2e{k) = H2(6), 



fc>i 



Note that Efc>iH2(^fc) e (5)_i,^(M) for every i < -1/2. 

We conclude the example with an observation that, although D^(iy) is not defined, 
Cy^^(W) is. If fact, f l3.16p implies that 

C^(W) = w, 

which is consistent with (12.141) and f l3.17p . 

If either u or w is a finite linear combination of ^a, then Du(y) is defined. The 
following proposition gives two more sufficient conditions for Du(f ) to be defined. 

Proposition 3.7. (1) Assume that there exist weights Tq,, ol ^ J such that 

J2 S'^'^cl^'clll < oo and ^ rlWu^Wl < oo. (3.18) 

If 

sup — — := bd < oo (3.19) 

for every a ^ J , then D„(f ) is well-defined and 

II (D„(.)) J|^«^ < 2l-l hi Y: 2\%'\\v,rx E -l\M\u ■ (3-20) 

f3ej I3ej 

(2) Assume that there exist weights r^, ex. E J such that 

J2 rilKfx < oo and 5^ 2l"lr^l««||^ < oo. (3.21) 

aeJ ctej 

If 

rp 
sup := Cq < oo (3.22) 

for every a G JT", then D„(f) is well-defined and 

II {^u{v))Sx^u < 2l"l cl J2 -l\Mx E ^'%'\M\l ■ (3.23) 

I3£J I3ej 

Proof. Using 

fc>0 ^ ^ 

we conclude that (^) < 2" for all k > and therefore 

;)^n(-),2H (3.4, 
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for 8l\\ j3 E J . Therefore, 



^u{v)) 



OL I I X®W 



^ V ( " /3 ) ^°+^ ® ^^ 



x<m (3.25) 



and the result follows by the Cauchy-Schwartz inequality. D 

Remark 3.8. (a) If r^ = q" for some sequence q, then both fl3.19p and ( I3.22p 
hold, (b) More information about the structure of u and/or v can lead to weaker 
sufficient conditions. For example, if {ua^upju = for ck 7^ /3, and ||Ma||i/ < 1? then 

|2 



(D«(^0)c,L®w < 00 if and only if 



pej 



which is a generalization of (12. 6p . Similarly, if (n^, n^);^ = for o; 7^ /3, then (£„(!;)) 
exists for all a. E J and 



(^"'"')«^S6) 



^'/3 WUoc^fiWli- 



The reader is encouraged to verify that 

(1) If M = W , with Me(fc) = Ufc and Mq = otherwise, then (I3.14p . (13.151) . and 
fIXTBD become, respectively, ([22D, (ESD, and fl^TT]) . 

(2) The operators 5^^ and D^^ are the creation and annihilation operators from 
quantum physics [2j: 

Dffc(^a) = A/«fc^a-e(fc), <5efc(^a) = V«fc + 1 {a+e(fc). (3.26) 

More generally. 



De.(e«) = w (^) e.-/3, '5,,(e.) = \i{^Y) e"+/3' A^(e.) = (^) e.. (3.27) 

(3) If 

vEL2i¥-X), fEL2{¥-X®U), D„(t;)GL2(F;X®W), (5„(/) G L2(F; X), (3.28) 

then a simple rearrangement of terms shows that the following analogue of 
(I23|) holds: 

E(D„(tO, /)^^^ = E(i;, 5„(/))^. (3.29) 

For example. 



D«(^7) = Xl y f j W7-"'^c 
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and, if we assume that u and / are such that <5„(/) G L2{¥; X), then 

aej V ^'^^ 

(4) With the notation Hq = VoJ.^a, 

and 

5h„(H^) = H„+/3. (3.30) 

To conclude the section, we use ( I3.30p to estabhsh a connection between the Sko- 
rokhod integral <5 and the Wick product. 

Definition 3.9. Let f be a generalized X-valued random element and rj, a general- 
ized real-valued random element. The Wick product f orj is a generalized X-valued 
random element defined by 




foV=2^(2^\l{a]fc.-f3V(3] e«. (3.31) 

a£j \l3ej 

The definition implies that / o ?/ = 77 o /, 

^a0^l3 = \ i 1 ^a+/3, Hq, O H^ = Ea+f3- (3.32) 

In other words, f l3.3ip extends by linearity relation f l3.30p to generalized random 
elements. Comparing f l3.3ip and f l3.15p . we get the connection between the Wick 
product and the Skorokhod integral. 

Theorem 3.10. If f is a generalized X-valued random element and rj, a general- 
ized real-valued random element, then Sn{f) = f orj. In particular, if rj and 9 are 
generalized real-valued random elements, then 

d,{e) = de{v) = voo. 

The original definition of Wick product [2D] is not related to the Skorokhod integral, 
and is it remarkable that the two coincide in some situations. The important fea- 
ture of fl3.3ip is the presence of point-wise multiplication, which does not admit a 
straightforward extension to general spaces. 

A natural definition of the multiple Wiener-Ito integral in the one- dimensional setting, 
that is, with respect to a single standard Gaussian random variable ^, is as follows. 
With only scalar as possible integrands, set 

/„(1)=H„(0. 
As expected, 

0^(4(1)) = n/„_i(l), 6^{In{l)) = In+lil). 

This is consistent with the general definition as long as the Wick product is used 
throughout: /n(l) = C,^"'- An interested reader can easily extend this construction to 
finitely many iid standard Gaussian random variables. 
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Definition 13.91 and Theorem 13.101 raise the following questions: 

(1) Is it possible to extend the operation o by replacing the point-wise product 
on the right-hand side of (13.311) with something else and still preserve the 
connection with the operator 67 Clearly, simply setting f ou = du{f) is not 
acceptable, as we expect the o operation to be fully symmetric. 

(2) Under what conditions will the operator v ^-^ uov be (a Hilbert space) adjoint 
or (a topological space) dual of D„? 

(3) What is the most general construction of the multiple Wiener-Ito integral? 

We will not address these questions in this paper and leave them for future investi- 
gation (see references [151 [16] fo^^ some particular cases). 

4. Elements of Malliavin Calculus on special spaces 

The objectives of this section are 

• to establish results of the type 

\\Aaiv)\\a<C{\\u\\b) \\v\\c, 

where || ■ ||j, i = a,b,c are norms in the suitable sequence or Kondratiev spaces, 
the function C is independent of v, and A is one of the operators D, d, C 

• to look closer at D and 5 as adjoints of each other when fl3.28p does not hold. 

To simplify the notations, we will write L2JX) for L2n(F;X). 

We start with the "path of the least resistance" approach and see what one can 
obtain with a straightforward application of the Cauchy-Schwartz inequality. The 
first collection of results is for the sequence spaces. 

Theorem 4.1. Let q = {qk, k > 1} be a sequence such that qk > 1 for all k and 
J2k>i Vo'fc < o*^- Denote by y^q the sequence {\/2qk, k > 1}. 

(a) Ifue L^liU) and v e L^^^^^X), then Bu{v) G L2{¥;X ®U) and 

(b) Ifue L-,l{U), f e L:^1{X®U), and Efc>i2V?fc < 00, then S^if) E L-^^^{X) 
and 

In particular, if u E L2^^{hl) and f G L^q(X (g) U), then 6u{f) G L2^^{X). 

(c) Ifue L^^iJA), V G ^2,v/2q(^). and Y.^^^ 2^/ql < 00, then Cu{v) G L-^^^{X) and 

(\ 1/2 / \ 1/2 
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Proof, (a) By (Km with r„ = 6^, = q^", 

The result then follows from (13.101) . 

(b) By (13.151) . fl3.24p . and the Cauchy-Schwartz inequality, 

/3 f3<cx 

and the result follows. 

(c) This follows by combining the results of (a) and (b). D 

Analysis of the proof shows that alternative results are possible by avoiding inequality 
(13.24^ : see Theorem 14.31 below. The next collection of results, this time for the 
Kondratiev spaces, is again in the spirit of the "path of the least resistance." 

Theorem 4.2. (a) If u e (5)„i,_^(W) and v G {S)i,i{X) for some i e R, then 
T>u{v) e {S)i^i_p{X ®U) for dip > 1/2, and 



1 ■'^' 

-l--l-r^72A 



|D„,M|iyc' _ ,.^,.^ < I I I Z ^7^^ I lk||,S^ , „fm \\V\ 



(5)i,,_p(X®W) ^ Mll_(2fc)-2P I ll"ll('5)-i,-HW) IK'll{5)i,,(X). 



^fc>l 



(b) If u e (5)_i,£(W) and f e (5)_i,^(X ® U) for some i e M, then d^if) e 
(iS)_i^^_p(X) for every p > 1/2, and 



\aej J 



1/2 

lkl|(5)_i,f(W) \\J \\{S)^ij{X®U)- 



In particular, if u G S^i{U) and f G iS„i(X <^U), then <5„(/) G 5_i(X). 

(c) If u & {S)-i-i{U) and v G (iS)i^+p(X) for some C. E M. and p > 1/2, then 
Cuif) e (5)_i,,_p(X) 



l/:.MII.... , „ „... < I Ut^^) (^'£{2Nr^-' 



2 



l"ll(-S)_i,_f(W) ll^ll(5)i,^(X®W)- 



Proof. To simplify the notations, we write r^ = (2N)^". 
(a) By dan. 

To get the result, use triangle inequality, followed by the Cauchy-Schwartz inequality 
and (13:91) . 
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(b) By fl3.15p and the Cauchy-Schwartz inequality, 

and the result follows. 

(c) This follows by combining the results of (a) and (b), because 

(5)i,,(X) C (5)_i,,(X). 



D 



Let us now discuss the duality relation between 5u and Vu- Recall that fl3.29p is 
just a consequence of the definitions, once the terms in the corresponding sums are 
rearranged, as long as the sums converge. Condition fl3.28p is one way to ensure the 
convergence, but is not the only possibility: one can also use duality relations between 
various weighted chaos spaces. 

In particular, duality relation (13. 5p and Theorem 14.21 lead to the following version of 
(IX^ : if, for some £ G M and p > 1/2, we have u G (5)„i,_£_p(W), v G {S\t+p{X), 
and / G (5)_i,^(X ® W), then 

(5„(/),t;)i,,+p=(/,D„(t;))i,,. (4.1) 

To derive a similar result in the sequence spaces, we need a different version of The- 
orem U]T1 

Theorem 4.3. Let p, q, and r he sequences of positive numbers such that 

111 
- + - = -, A; >1. (4.2) 

PI Tk n 

(a) Ifue L2,p{U) and f G -^2,q(X (S)U), then du{f) G -^2,r(-^) and 

\\^uif)\\L2,AX) < \\u\\l2,p{U) ||/||L2,q(X®W)- 

(b) In addition to (14. 2 p assume that 

Ei<oc. (4.3) 

k>iPk 



Define 



\ 1/2 



Ifue -^2,p(W) and V G L^^iX), then Bu{v) G L^^^iX ®U) and 

l|D«(tOllL-J(x®w) < C'||n|U,_p(w) llt'IL-j(x)- 
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Proof, (a) By (l3J[5ll . 



lum 



2 



sE 

-rej 



E 

ci+/3=7 

E 

c»+/3=7 



7 



q; 



{U,Uf3) 



U 



r27 



X 



|(/o,W/3)w|t"r^ 



X 



:i + c-')°r'" = q2", (1 + c)°r2" = p'". 



Define the sequence c = {c^, /;; > 1} by c^ = pl/ql, so that 
Then 

II t- / n\ 11") , X ^ I \ ^ I I ^ \ rv/9 _/->./9ll ^^ 



(4.4) 



i''.(/)iii,,.,x, < E f E \/(l)'°'''"'""ii/" 



a^/? 



76J' \q;+|3=7 

By the Cauchy-Schwartz inequahty and (13.121) . 



I^.(/)II 



L2,r{^) 



< 



E E 



E '-''Wfc^WmxWMgv'-t 



2/3 



\a+l3=-f 



l|2 II ||2 „2a 2/3 



76^" V \a+/3=7 

7eJ'a+/3=7 

E wf^Wu^x (1 + c-^)v-) (e \Ml (1 + ov^) 

^aej / \f3ej J 



\aej J \i3ej 



2/3 



2 Il7/I|2 

L2,p{X(S>U) II"IIL2,p(W) 



(b) By dSH, 



|D«Miii~j(x«,,) = E 



aej' 



E 

/3GJ' 



Ct + 

(3 



Va+I3 ® ^^/3 



-2a. 



X(g)U 



a + /3\ , _2a 

^ ) h'»+(3\\x \M\u I q 



^E E 

aej \l3ej 
Define the sequence c = {c^, fc > 1} by c^ = r^/pl < 1, so that 



(c-^ - l)°q2" = p2°^ (1 _ c)"q2" = ,:2«_ 



(4.5) 



Then 



|D4^)lli-j(x«^,) 



< 



E E 

aej \i3ej 



« + /3') c/^/2q-(«+/^)c-/3/2||^^^^||^q/3||^^||^ 



V ,«l|2 r-/3n-2(a+/3) 2/3 ii ||2 
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By the Cauchy-Schwartz inequality and (13.101) . 

iiD.('')iii-;„«„, < Z (e ("a")'") f5:'!-**'"'-''llwR".^''ll".ie 

a6j V /36J 

/3GJ \ \ae:7 / / 

\/3ej7 / \a6j 

-r^ii ||2 II ||2 

where the last equality follows from (14. 5p . Note also that 

En , ||2 -2(a+/3) ^ n ,||2 

and the equality holds if and only if /3 = (0) . D 

Together with duality relation (13.41) . Theorem 14.31 leads to the following version of 
(I3.28P : if n G L2,p(W), / G -^2,c|(-^ ®U), and v E L2^{X), if the sequences p,q,t are 
related by (g^D, and if (H^D holds, then 

{Su{f),v),= {f,B^{v)),. (4.6) 

Here is a general procedure to construct sequences p,q,r satisfying (14. 2 p and (14. 3p . 
Start with an arbitrary sequence of positive numbers p and a sequence c such that 
< Cfc < 1 and J2k>i ^k < oo. Then set r^ = C}^p\ and q\ = pi/{c^^ ~ !)• If the space 
U is ra-dimensional, then condition (14.30 is not necessary because sequences p, q, r are 
finite. 

Theorem 4.4. Let p, q, anc? r he sequences of positive numbers such that 

p^g^ > 1, A; > 1, and E ^^ < ^- (4-8) 

fc>i ^fc^fc 

If u E L2,p(W) anc? t> G L2,q(X), i/ien Cuiy) G L2,r(-^) o^t? 

(\ 1/2 
n;#h) ll«llL,p(t.) II^IU..m. (4.9) 
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Proof. It follows from f l3.16p that 



l(£.W)Jli<(E/( 



) (/^jll^^+^ll^ hc.-f3\\u \\U^\\U 



Let f) = {hk, k > 1} he a sequence of positive numbers such that that 

/lyfc < 1, A; > 1, 2. ^k < oo. 



k 

Define 



^- ni^ 



1/2 

- hk 



Then 

1/2 / \ 1/2 






-711 ||2 

■ \/ \ T / \ \ T / / \ 

7 

1/2 



Cft, I t: 7T77 1 I > fl '^ll'«fl+7llxll'"7L 

Next, take another sequence tti = {wk, A; > 1} of positive numbers and define the 
sequence c = {ck, k > 1} hj 

Ck = j^. (4.10) 

1 -hk 

Then 



As a result, 

^,2(/3+7)(l + c)^+'^rt)-(^+^)||t;;3+7l|^ 
/3 

Then (gS]) holds if 

rl{l + Ck)= ..7 ^, =pI, '^^^^±^ = ql (4.11) 

rfc(l + Ck)hk Wk 

The three equations in (14. lip imply 

/1 \ Pk Pk 1 

(,1 + Ck) = —n, Wk = ^, Ilk = ~2~2' 

'"fc Qk PkQk 
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and then ( 14. 7p follows from f l4.10p . Note that a particular case of (14 .Zp is qk = I/t/c, 
p^^ + 1 = r^^, which is consistent with Theorem 14.31 if we require the range of D^ to 
be in the domain of 6u- 



D 



Example 4.5. Let W = X = M. Then a = n e {0, 1, 2, . . .}, 



n>0 n>0 n>0 

To begin, take u = ^. Then 

n>l n>0 

n>l 

Next, let us illustrate the results of Theorems 14.31 and 14.41 Let p, q, r be positive real 

numbers such that 

1 1 _ 1 

p q r ' 



for example, p = q = 1, r = 1/2. By Theorem 14. 3[ if 

Y.p'^uIkoo, $^;f <oo. 



n>0 n>0 



then 



and if 



then 



n>l ^ 

Xp"n2 <oo, J]g"/2<oo, 

n>0 n>0 

n>l 

If p, q, r are positive real numbers such that 

l_i)(,_i 

r pj \ p 
(for example, p=l, q = 2, r = 1/2) and 

n>0 n>0 

then, by Theorem 14.41 

J2r-{C^{v))l <oo. 

n>0 
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